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Abstract 

> 

\ We describe the current search for confinement mechanisms in lattice QCD. We 

report on a recent derivation of a lattice Ehrenfest-Maxwell relation for the Abelian 
projection of SU(2) lattice gauge theory. This gives a precise lattice definition of field 
qq | strength and electric current due to static sources, charged dynamical fields, gauge 

fixing and ghosts. In the maximal Abelian gauge the electric charge is anti-screened 
analogously to the non-Abelian charge. 
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1 Introduction 



Quenched lattice QCD calculations of the static quark anti-quark potential 
have firmly established a linearly rising behavior over all distances obtain- 
able in state-of-the-art simulations^). However this situation satisfies no one, 
least of all a Dick Slansky. A 'black box' calculation at limited quark sepa- 
ration gives supporting evidence that QCD confines quarks, but it offers no 
explanation nor reveals a principle governing the phenomenon. 

Lattice QCD is more than just an algorithm to calculate quantities at strong 
coupling. The lattice is a regulator for QCD, parametrized by the lattice spac- 
ing a. Gauge symmetry is preserved at all costs. This regulation scheme is 
perhaps the only one that gives a completely self consistent cutoff model in its 



1 Talk given at the Symposium in Honor of Dick Slansky, Los Alamos, May 20-21 
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own right. The dynamical variables are elements of the symmetry group rather 
than the Lie algebra. As a consequence, many of the topological features that 
are 'likely suspects' in the physics of confinement have natural definitions on 
the lattice. These include conserved U(l), Z(N), SU(N)/Z(N) monopole loops, 
Dirac sheets, Z(N) and SU(N) vortex sheets and other features. 

Lattice work(|2]; |3|; |j; ^ is coming into its own in sorting through some of the 
seminal ideas on confinement (0; |8|; |9|; [10]; [11]; ]12|) proposed in the '70's. I venture 
to say that Dick Slansky could easily have been drawn into the fray because 
it has all the elements that typically captured his imagination. The central 
questions are fundamental. They involve interesting issues in group theory, 
topology, and duality. Candidate mechanisms are proliferting and fundamental 
questions remain unsettled^ 

Confinement is a consequence of a disordered state characterized by the ex- 
pectation value of a Wilson loop suppressed to an area law rather than a 
perimeter law for an ordered system. Contributions from topological objects 
e.g. monopoles and vortices can accomplish this. On the lattice these objects 
are not singular. They occur abundantly in SU(N) lattice theories. It is only 
in the continuum limit of zero lattice spacing where these approach singular 
structures. In lattice simulations, one can tamper with these objects and see 
if suppression is accompanied by the loss of an area law. In this way one has 
a laboratory to study candidate disordering mechanisms. 

It is not realistic to try to give a proper review of the many competing views in 
this brief article. Further, the studies are possibly diverging more rapidly than 
at any point in the past considering the variety of papers on the subject over 
the past year. However I would like to touch on them and further to argue that 
all the studies are describing genuine properties of QCD, seen through the eyes 
of subsets of the full dynamical variables. The controversies have to do with 
which scenario will lead to the most compelling explanation of confinement. 

Although these lattice studies reveal relevant confinement physics, the goal is 
not completely clear. Most would agree however that if we knew the dual form 
of QCD it would give a definitive description of the physics of confinement. 

Consider the example of superconductivity. By identifying the carriers of the 
persistent current one can discover an instability in the normal vacuum. The 
Ginzburg-Landau effective theory describes the consequence of the instabil- 
ity which of course elucidates the fundamental principle underlying the phe- 
nomenon which is the spontaneous breaking of the U(l) gauge symmetry ([13|). 
We might imagine that an understanding of the topological structures might in 

2 At Lattice '98, the XVI International Symposium on Lattice Field Theory in 
Boulder, there were 52 contributed papers on color confinement, up from a dozen 
or so in the early '90 's. 
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the same way lead to a discovery of a spontaneous gauge symmetry breaking. 
The problem is that a definitive scenario still seems quite distant. 



2 Lattice confinement studies - a brief survey 



In this section I wish to call attention to some of the approaches in lattice 
studies of confinement that have been active over the past 12 months. This is 
principally to give a flavor of the work and a list of recent references, and rely 
in the references therin for a more complete bibliography. 

Our most recent work, which we describe in Sec. 5 is based on the Abelian 
projection. 

We will restrict our attention to SU(2) theories throughout this article, be- 
lieving that the essential issues will be revealed in this simplest case. 



2.1 Abelian projection, Abelian dominance 



In the Georgi-Glashow model with a gauge field coupled to the adjoint Higgs 
one can define a gauge invariant Abelian field strength (|14]). This is the con- 
struction used to identify the magnetic field of a 'tHooft-Polyakov monopole. 
Further, the Higgs field can be used to define an Abelian projected theory. 
This is accomplished by gauge transforming the Higgs field into the 3 direc- 
tion. The resulting partially gauge fixed theory is still invariant under U(l) 
gauge transformations - rotations about the 3 axis. 

Pure gauge SU(2) has no adjoint Higgs field and so there is no straightfor- 
ward way to define Abelian variables. It is still possible to define an Abelian 
Projection^; [Hj; 0). 



One can find a "collective" Higgs field transforming under the adjoint repre- 
sentation which we denote as an "adjoint field". Consider an arbitrary Wilson 
line starting and ending at a particular site. This can be parametrized by 
W(x) = cos x + 20 ■ cr sin \ where (p is a normalized adjoint field. The Wilson 
line could be (i) an open (no trace) Poliakov loop, (ii) an open plaquette in 
e.g. the (1, 2) plane, or (iii) an arbitrary sum of such lines then normalized to 
construct an SU(2) element. One can generalize to (iv) define the adjoint field 
self consistently by requiring that the adjoint field at one site be equal to the 
normalized sum of the eight neighboring adjoint fields after parallel transport 
to that site. 

In these four examples, one can further fix the gauge by rotating the adjoint 



3 



field into the third direction. This then fixes (i) to the Polyakov gauge, (ii) to 
a plaquette gauge (iv) e.g. a clover gauge, and (iv) in the maximal Abelian 
gauge. The "collective" adjoint field itself is gauge covariant, it does not neces- 
sitate a gauge fixing. Different adjoint fields correspond to different Abelian 
projections. 

The maximal Abelian gauge appears to be the most promising choice. The 
condition, equivalent to the one given just above(|T7|), maximizes X)(t/| + U\ — 
C/f — Uf), where U = U4 + icr ■ U , U\ + U 2 = 1. The maximization brings the 
links as close as possible to the values Ui = U% = and JJ\ + C/| = 1, leaving 
a U(l) invariance. 



A new and interesting alternative has been proposed recently by Van der 
Sijs(|T9|). This gauge fixing is based on the lowest eigenstate of the covariant 
Laplacian operator: i.e. sum the adjoint fields at the eight neighboring sites - 
parallel transported to a given site - and subtract 8x the field at the site. This 
defines a gauge with no lattice Gribov copies. And it identifies the positions 
of monopoles as singularities (zeros) of the adjoint field. 

'Abelian dominance' asserts further that operators built out of the U(l) links 
defined by the 'Abelian projection' will dominate the calculation of string ten- 
sion and other observables. There are a large number of numerical tests which 
support Abelian dominance. However the successful tests are for the most part 
restricted to (a) certain specific quantites e.g. string tension in the fundamen- 
tal representation, and (b) for the gauge choice of the maximal Abelian gauge. 
In other gauges, e.g. plaquette gauge, and/or for other quantities such as string 
tension in other representations Abelian dominance is not well supported and 
in some cases strongly violated. 



On the other hand in a series of papers Di Giacomo et. al.( f20! pT| p?2j ; |23| ; [24]) 
have studied the finite temperature confining phase transition and shown that 
the disorder parameter is very insensitive to the choice of Abelian projection. 
This is the foundation of the dual superconductor scenario. This is not in con- 
flict with the above results because it does not require the Abelian dominance 
ansatz. 



Using the Abelian projection and Abelian dominance, the problem is reduced 
to a U(l) gauge theory albeit with a complicated action coming from the ef- 
fects of gauge fixing. As in a pure U(l) gauge theory, Dirac monopoles can be 
identified as the charge carriers of the persistent currents in the dual super- 
conductor scenario^; |IBJ g§ ^). 

Bakker, Chernodub and Polikarpov(j30l) have shown that Abelian monopole 
currents, defined in the maximal Abelian gauge are physical objects: there is 
a strong local correlation between monopoles and enhancements in the SU(2) 
gauge invariant action. 
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Other recent results of Polikarpov et al(^TJ [32]; [33]) include: Abelian monopoles 
in the maximal Abelian gauge are dyons. There are strong correlations between 
magnetic charge, electric charge and topological charge density. These same 
connections have been reported independently by Ilgenfritz et al(]34]). 



In order to understand the disordering effect of monopoles on large Wilson 
loops, Hart and Teper(p5|) have studied the clustering properties of monopole 
loops and found two classes of clusters: A single cluster permeates the whole 
lattice volume, the remaining are small localized clusters. 



In a series of papers, Ichie and Suganuma(PB| |5T| ; have sought a deeper 

understanding of Abelian dominance. They look at the residual degrees of 
freedom after Abelian projection, i.e. the coset fields and argue that in a 
random variable approximation that Abelian dominance is exact. 



Ogilvie(|40|) has argued in the context of the Abelian projection that gauge 
fixing is in principle unnecessary, that results are the same whether the gauge 
is fixed or not. This is at odds with many simulations, leaving much to be 
sorted out. 

Grady([4ll |42]) has given evidence that casts doubt on whether Abelian monopoles 
confinement mechanism carries over to the continuum limit. 

A number of authors, including our group have reported that a well defined 
electric vortex forms between static sources(]33|; |T7]; fP|). The vortex is very 
well described by a dual Ginzburg-Landau effective theory. In effect we are 
able to show that there is a local relation between the electric flux and the 
curl of the monopole current, defined by Abelian projection of full SU(2) 
field configurations that gives a damping of fields as they penetrate the dual 
superconductor. This directly accounts for confinement. In Sec. 5 we derive 
a precise definition of the electric field which tightens up the definition of 
vortices 



2.2 Vortices in Z(N) x SU(N)/Z(N) formulation 



In the early 80's there was an effort by Mack and Petkova(|46|), Yaffc(^7|), 
Tomboulis (f£8|), Yoneya(|49|), Cornwall (|50|), and Halliday and Schwimmer(|5T|), 
to understand confinement in terms of SU(N) vortices. These are singular 
gauge configurations characterized by their topological properties. They are 
multiple valued in SU(N) but single valued in SU(N)/Z(N). 

There are more recent developments by Tomboulis and Kovacs (p2| [53]; 0) [[]. 

3 A thorough review of the topological issues in this approach can be found in 
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Their study is based on a decomposition of the SU(2) partition function into 
variables defined over the center Z(2) and variables defined over the quotient 
group SU(2)/Z(2) (= SO(3))f|). Whereas the SU(2) manifold has a trivial 
topology, SO(3) does not: tti(SO(3)) = Z(2). 

To accomplish this reformulation, one appends new Z(2) valued variables a p , 
which live on plaquettes, to the usual set of SU(2) links. This allows one to 
write the standard SU(2) Wilson action partition function in terms of these 
Z(2) variables, a p and link variables, {Ui}, defined over the manifold of SO(3) 
rather than SU(2). The latter distinction is manifest in the reformulated par- 
tition function in that the SU(2) links {U{\ occur only in Z(2) invariant com- 
binations i.e. invariant under the transformation Ui — > ±[7;. 

Consider a single (1,2) plaquette with a p = — 1 , all others = +1. A Z(2) 
monopole occurs in each of the two (1,2,3) cubes adjacent to the face of the 
negative plaquette and similarly in the two (1,2,4) cubes. 

On the dual lattice, each cube becomes a dual link orthogonal to the cube and 
this forms a conserved Z(2) monopole current loop (the smallest such loop) on 
the dual lattice. The negative plaquettes become dual plaquettes on the dual 
lattice. They form a surface bounded by the monopole loop. This surface is 
called a "thin vortex" sheet. 

One can build larger monopole current loops and vortices by laying out stacks 
of negative a plaquettes. 

There are also the usual plaquettes built from the trace of 4 SU(2) link vari- 
ables. Again consider a configuration in which all such plaquettes are positive 
except for a single (1,2) plaquette which is negative. An SO(3) monopole oc- 
curs in each of the two (1,2,3) cubes adjacent to the negative plaquette and 
similarly in the two (1,2,4) cubes. This is an SO(3) object because every link 
involved in its construction occurs quadratically and hence is invariant under 
a Z(2) transformation of each link. 

Analogously there is an SO (3) monopole current loop on the dual lattice. The 
surface bounded by the loop is a a Dirac sheet which differs from the "thin 
vortex" sheet in that it can be moved arbitrarily by multiplying the SU(2) 
links by elements of the center ± without costing any action. 

The equivalence of the two forms of the partition function requires that every 
Z(2) monopole be coincident with an SO (3) monopole and vice versa. Hence 
we can visualize the vortex structure as two currents on coincident loops, 
bounded by two surfaces. This is called a 'hybrid' vortex. We can also shrink 
the monopole loop to zero, leaving either a pure "thin vortex" sheet, or in the 

Ref.(||) 
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other case a pure Dirac sheet. 



The SO (3) monopoles also have associated non-Abelian (non-integer valued) 
flux which costs action. These configurations are denoted "thick vortices", The 
Dirac sheet is the return flux required by the topological group iri(SO(3)) = 
Z{2). 

In this picture, it has been shown that thin vortices are exponentially sup- 
pressed at weak coupling and can not account for confinement. However SO(3) 
thick vortices can occur at weak coupling in which all terms in the action 
\trU p \ ~ 1. In recent papers Kovacs and Tomboulis(52; 53) have shown that 



the static quark potential can be reproduced by contributions from SO (3) thick 
vortices linking the Wilson loop in SU(2); and in SU(3); and that exclusion of 
vortices results in a perimeter law. 



See Gavai and Mathur(E3) for a study of Z(2) monopoles and the deconfine- 



ment phase transition. Also see Grady (|56|) for variation on the SO (3) - Z(2) 
monopole construction. 



2.3 Center Projection Vortices 



There are other closely related approaches that also focus on the center of 
the group as the key to understanding confinement. In a number of recent 
papers, Greensite, Del Debbio, Faber and Olejnik(^) have introduced the 
'Center projection' as a way of identifying Z(2) vortices for SU(2). which they 
denote as "Projection vortices". They differ from the thick and thin vortices of 
Tomboulis. They are defined in the the maximal center gauge. In this gauge, 
one maximizes J2U%. Hence U4 will be as close as possible to ±1. The Z(2) 
links are then defined Z = sign(Ui). This gives a Z(2) gauge theory that has 
"P vortex" excitations which are somewhat analogous to the Tomboulis 'thin' 
vortices. They also define "center vortices" which can be much larger than one 
lattice spacing and are analogous to the Tomboulis "thick vortices". 



Greensite et. al. fl58[; |59t |60|) have reported a growing list of encouraging results. 
Among the results are: P vortices are well correlated with SU(N) vortices; No P 
vortices no confinement; P vortices account for the full string tension; P vortex 
density scales; Center vortices thicken as the lattice cools; P-vortex locations 
are correlated among Gribov Copies; Preliminary successful generalizations to 
SU(3); and center vortices are compatible with Casimir scaling. 



See also Langfeld et alflBTt |62] ) and Stephenson (|63|) for further results on center 
vortices. 

A variation in the center projection procedure is to first fix to the maximal 
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Abelian gauge, and then maximize J2 cos X where x is the U(l) link angle. 
This is denoted the indirect maximal center gauge, as opposed to the above 
procedure denoted the direct maximal center gauge. In this gauge, a sheet 
consisting of monopole loops alternating with anti-monopole loops coincides 
with the P vortex sheet (|58|; |60| ) indicating a possible overlap between this and 
other scenarios. 



2.4 Dual variables 

One expects that the disorder regime in QCD would be described by an ordered 
regime in dual QCD. A definitive form of dual QCD would probably take us 
a long way toward an understanding of confinement. 

We would like to call attention to recent lattice work on dual variables by 
Cheluvar a j a (|6^) , and continuum work by Majumdar and Sharatchandra(|65*D , 
Sharatchandra et al(^), Mathur(|67|), Faddeev and Niemi(|68|), and Chan and 
Tsun(||). 

The work by Majumdar and Sharatchandra supports the dual QCD ansatz by 
Baker, Ball and ZachariasenflTOD . 

Seiberg and WittenflTTD exploited duality to establish confinement for super- 
symmetric QCD. 



2. 5 Instantons 

An important aspect of confinement studies is to identify what objects can 
not account for confinement. Instantons were likely suspects at one time. 
Smoothing techniques are very important in the identification of instantons. 
Unlike monopoles, small instantons can fall through the holes in the lattice 
and further can be swamped by short distance fluctuations. New smoothing 
techniques which overcome these difficulties, denoted 'renormalization group 
mapping', have been applied to configurations by DeGrand, Hasenfratz and 
Kovacs ([72|) to elucidate the role of instantons in the QCD vacuum. The come 
to the strong conclusion that instantons alone do not confine. 

See also other recent works by Narayan and Neuberger(|73|); Narayan and 
Vranas([F4]); Ph. de Forcrand, M. Garcia Perez and I.-O. Stamatescu(|75|); and 
B. Alles, M. D'Elia and A. Di Giacomo(|76]). 
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2. 6 Interconnections 



Even though the possible explanations for confinement seem to be diverging 
at this point, I would like to re-emphasize that different descriptions can in 
some cases be describing the same underlying physics: 

Abelian monopoles in the maximal Abelian gauge correlate with 

• SU(2) enhancements in the gauge invariant action fl30|), 

• gauge invariant topological charge (PT|; P^ ; 

• P vortices in the maximal center gauge (|59|), 

• Instantons([77|; |7|). 

Therefore the monopoles can not be thought of as merely gauge artifacts. 



3 U(l) gauge theories and superconductivity 



The onset of superconductivity is governed by the spontaneous breaking of the 
U(l) gauge symmetry via a non-zero vacuum expectation value of a charged 
field fljjD- An immediate consequence of this is the generation of a photon 
mass and, for type II superconductors, the formation of magnetic vortices 
which confine magnetic flux to narrow tubes ([751) as revealed by the Ginzburg- 
Landau effective theory. Lattice studies of dual superconductivity in SU(N) 
gauge theories seek to exploit this connection in establishing the underlying 
principle governing color confinement. 

In U(l) lattice pure gauge theory (no Higgs field), this same connection is seen 
to be present, not in the defining variables, but rather in the dual variables. 
More specifically: 

(1) A field with non-zero magnetic monopole charge, $, has been construc- 



ted(80D. It is a composite 4- form living on hypercubes constructed from 
gauge fields. There are also monopole current 3-forms. On the dual lattice 
this monopole operator is a 0-form living on dual sites. The monopole 
currents are 1-forms living on dual links. These currents either form 
closed loops or end at monopole operators. The monopole operator has 
a non-zero vacuum expectation value in the dual superconducting phase, 
($) 7^ 0, thereby signaling the spontaneous breaking of the U(l) gauge 
symmetry. 

(2) Dual Abrikosov vortices have been seen in simulations (E7I; BJ). They 
are identified by the signature relationship between the electric field 
and the curl of the monopole current in the transverse profile of the 
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vortex. The dual coherence length, ^ measures the characteristic dis- 
tance from a dual-normal-superconducting boundary over which the dual- 
superconductivity turns on. The dual London penetration length, A^ mea- 
sures the attenuation length of an external field penetrating the dual- 
superconductor. The dual photon mass ~ 1/A^ and the dual Higgs mass 

A signal (<&) 7^ without the consequent signal of a dual photon mass does 
not imply confinement. An observation of a dual photon mass, i.e. vortex 
formation, without (<&) 7^ does not reveal the underlying principle governing 
the phenomenon. 

3.1 Higgs effective theory 

The Higgs theory, treated as an effective theory, i.e. limited to classical so- 
lutions, and considered in the dual sense, provides a model for interpreting 
simulations of the pure gauge theory that can reveal these important connec- 
tions. Recalling the Higgs' current 

IP 

^ = -3 (<f>* ~ ieA ») <t>-4>{d» + ieA,) 0*) , (1) 

and spontaneous gauge symmetry breaking through a constraint Higgs poten- 
tial 



= ve ieuj{x \ v = constant, (2) 

leads to the London theory of a type II superconductor. Using Eqns.(l) and 
(2) we obtain 



J M = -eV (A, - d,cu) , (3) 
(df, J v - d u J^) + (d^J u - d u J^) = 0, 

V x J + ^-B = 0, 

A z 

where 



2 2 2 1 / A\ 

e v = m 7 = — . (4) 
Using Ampere's law V x B = J, we obtain 
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V 2 B = B/X 2 , 



identifying A as the London penetration depth. 

If the manifold is multiply connected, then the gauge term in Eqn.(3) can 
contribute, as long as eu(x) is periodic, with period 2n on paths that surround 
a hole. 



J (B+X 2 V x J) • nda = j(A+X 2 J) ■ dl, 

s c 

<j> Vw • dl, 



c 
2vr 

= N— = Ne m = $ m . 

e 

where N quanta of magnetic flux penetrates the hole in the manifold. In real 
superconductors, the hole is a consequence of the large magnetic field at the 
center which drives the material normal. 

A cylindrically symmetric vortex solution is given by 



B + A 2 V x J = $ m 5 2 (r ± )n z , 
(l-X 2 V 2 ± )B z (r ± ) = $ m 5 2 (r ± ), 

B z (r ± ) = ^K (rJX). 

The delta function core of this vortex is normal, i.e. no spontaneous symmetry 
breaking, and the exponential tail of the vortex is a penetration depth effect at 
the superconducting-normal boundary. The key point is that the modulus of 
the Higgs field must be independent of position to get these idealized vortices. 
For a "Mexican hat" Higgs potential, there is a coherence length setting the 
length scale from a normal-superconducting boundary over which the vacuum 
expectation value of the Higgs field changes from zero to its asymptotic value. 

On the lattice, the same phenomenon occurs. We can generate vortices from 
finite configurations. In the continuum these objects are singular. Since the 
lattice formulation is based on group elements, rather than the Lie algebra 
the periodic behavior of the compact manifold is manifest. This gives the 
2nN ambiguity in the group angle leading to N units of quantized flux. To see 
how this works, consider the lattice Higgs action 
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S = f3 Y, (l-cos9^(x)) 

X,fl>lS 

-KY{^{x)e i6 ^4>{x + n) + H.c.) + Y V m99S (\<p(x)\ 2 ), 

X,fl X 

where 9^ v {x) is the curl of the gauge field, 



and where (j){x) is the Higgs field and <fi(x + //) refers to the Higgs field at the 
neighboring site in the /i direction and A+ is the forward difference operator. 
The electric current is given by 



n 3 

—J e (x) = Im{ct)*{x)e id ^ x) (t){x + n)), 

where a is the lattice spacing. Let us choose a Higgs potential that constrains 
the Higgs field = 1. Then if 

sin[fl + 2iV7r] « 6, 

we find a relation between the field strength tensor and the curl of the current 



a 2 (A+Jfo) - A+J*(x)) _ 2nN 1 __ AT 1_ 
e^K a e cr 



where 



ea 2 ^ = sin[6> M (x) + 9 v (x + /i) - 6» M (x + i/) - ^(x)]. 

If AT = then this is a London relation which implies a Meissner effect. 
For iV ^ there are iV units of quantized flux penetrating that plaquette, 
indicating the presence of an Abrikosov vortex. 



3.2 Pure U(l) gauge theory 



In a pure U(l) lattice gauge simulation (without Higgs field), lattice averages 
over many configurations exhibit superconductivity in the dual variables. The 
superconducting current carriers are monopoles. They can be defined in a nat- 
ural way on the lattice using the DeGrand Toussaint (p6|) construction. They 
arise in non-singular configurations again because of the 2nN ambiguity in 
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group elements. These are the magnetic charge carriers for dual superconduc- 
tivity. For a review of the monopole construction and vortex operators in U(l) 
gauge theory see e.g. the 1995 Varenna Proceedings([79|) . 

As a brief summary, consider the unit 3- volume on the lattice at fixed X4. The 
link angle is compact, —71 < 6^ < ir. The plaquette angle is also compact, 
— 4-7T < 9 av < 4ir and defined 



ea 2 F^ u {x) = 6^(x) = &+6 v (x) - A+^(x), 

where a is the lattice spacing. This measures the electromagnetic flux through 
the face. Consider a configuration in which the absolute value of the link 
angles, 16^1, making up the cube are all small compared to 7r/4. Gauss' theo- 
rem applied to this cube then clearly gives zero total flux. Because of the 2n 
periodicity of the action we decompose the plaquette angle into two parts 



0^{x) = 9^ u (x) + 2rrn flu (x). (6) 

where — it < 6^ v < ir. If the four angles making up one of the six plaquette 
are adjusted so that e.g. 6^ v > tc then there is a discontinuous change in 
9 fiv by — 27r and a compensating change in n^. We can clearly choose the 
configuration that leaves the plaquette angles on all the other faces safely 
away from a discontinuity. We then define a Dirac string n^ u passing through 
this face (or better a Dirac sheet since the lattice is 4D). M „ measures the 
electromagnetic flux through the face. 

This construction gives the following definition of the magnetic monopole cur- 
rent. 



Cm 



(7) 



This lives on dual links on the dual lattice. Although Eqn.(6) is not gauge 
invariant, Eqn.(7) is. Further it is a conserved current, satisfying the conser- 
vation law A+J™(x) = 0. 

In simulations of a pure U(l) gauge theory we find that lattice averages give 
a relation similar to Eqn.(5), but in the dual variables 



A-J^(x) - A~J™(x)) 1 

(X CL 



where *F flv is dual of F^. This is the signal for the detection of dual vortices([79|). 
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4 Non-Abelian theory 



The link of these considerations to confinement in non-Abelian gauge the- 
ory is through the Abelian projection |TB|). One first fixes to a gauge 
while preserving U(l) gauge invariance. The non-Abelian gauge fields can 
be parametrized in terms of a U(l) gauge field and charged coset fields. The 
working hypothesis is that operators constructed from the U(l) gauge field 
alone, i.e. Abelian plaquettes, Abelian Wilson loops, Abelian Polyakov lines 
and monopole currents, will exhibit the correct large distance correlations rel- 
evant for confinement. 



In the continuum limit the maximal Abelian gauge condition is 



(d,±gAl(x)) A±(x) = 0. 



This is achieved on the lattice by a gauge configuration that maximizes R, 
where 



and where U^n) is the link starting a site n and extending in the /x direction. 



cos 



{<j>„{n))e i9 M sin^H^M \ 



sin(0 M (n))e- z ^(^) cos^n^e -4 '^ 71 ) 



After gauge fixing, the SU(2) link matrices may be decomposed in a 'left coset' 
form: 



cosf^^n)) sin(0 A1 (n))e l7M< ^ n ^ 



\ / e i0 M (n) q 



v -sin(0 M (n))e 



cos(<^(n)) ) 



\ 



Q e ~id^{n) 



Under a U(l) gauge transformation, {g(n) = exp [ia(n)o"3]}, 



0»(n) -> O^n) + a(n) - a(n + p) 7 M (n) -> 7^(71) + 2a(n) 



(9) 
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In other words, the left coset field derived from the link U^n) is a doubly 
charged matter field living on the site n and is invariant under U(l) gauge 
transformations at neighbouring sites. 

The c M = cos(0 M ) are real-valued fields which near the continuum ~ l + 0(a 2 ) 
where a is the lattice spacing. The off-diagonal = sin(0 A1 )e l7 ' J become the 
charged coset fields gaW^x), and 9^ the photon field gaA 3 ^(x). [The SU(2) 
coupling (3 = p- in 3+1 dimensions.] 

The static potential constructed from Abelian links gives as definitive a signal 
of confinement as the gauge invariant static potential as found by Suzuki 
et.al.(|n§ f§), Stack et.al.(^) and Bali et. al.(g). Bali et.al. find the Abelian 
string tension calculation gives 0.92(4) times the full string tension for (3 = 
2.5115. Whether this approaches 1.0 in the continuum limit remains to be 
seen. 



Equation (4) gives the connection between the non-zero vacuum expectation 
value of an order parameter and the photon mass or equivalently the London 
penetration depth. 

This connection between order parameter and penetration depth is the key 
to connecting spontaneous symmetry breaking to vortex formation and hence 
confinement. Dual superconductivity studies seek to establish a connection, of 
course calculated from the original variables. 



Di Giacomo et.al. 



21; 22; 23 



have reported extensive studies of the 



behavior of the order parameter for the dual theory (denoted disorder param- 
eter) at the confining transition. See also Polikarpov et. al.(^). Our group 
and others have reported dual vortex formulation between sources allowing 
the determination of the London penetration depth ([43]; O; 041). 



We can demonstrate a qualitative correspondence between these two indica- 
tors of dual superconductiviy in that they are both observed and both show 
the correct behavior on the two sides of the transition. However technical 
difficulties have eluded a direct comparison check of the dual form of Eqn.(4) 

Figure 1 shows a plot of p = 2^ In (//) as a function of (3 near the transition. 
The spike indicates step behavior in (/a) at the transition. 

Electric dual vortices between sources are well established(|43"l; |i~7| ; 44) in this 



gauge. The typical behavior is shown in Fig. 2. The London relation is seen 
in the confining case for transverse distances larger than about two lattice 
spacing. The dual coherence length ^ « 2, i.e. the onset of the violation of 
the London relation. The unconfined case is also shown where curl J is almost 
zero, and the electric field falls more gradually than in the confining case. 
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Bali et al(f44|) have done a large scale simulation of these vortices. Again fitting 
the electric field and the monopole current to the Ginzburg-Landau theory 
their results are shown in Fig 3 and 4. Their data is very well described by the 
two G-L parameters A = 1.84(24) and £ = 3.10(40) The ratio G-L parameter 
k = A/£ = 0.59(14). They find the total flux in the vortex, $ = 1.10(2) 
in units of the quantized flux. They concluded tentatively that type I dual 
superconductivity is indicated. 



5 Definition of electric field strength 

Central to finding the effective theory is the definition of the field strength op- 
erator in the Abelian projected theory, entering not only in the vortex profiles 
but also in the formula for the monopole operator. All definitions should be 
equivalent in the continuum limit, but use of the appropriate lattice expression 
should lead to a minimization of discretization errors. 

In a recent paper (|45|) we exploited lattice symmetries to derive such an oper- 
ator that satisfies Ehrenfest relations; Maxwell's equations for ensemble aver- 
ages irrespective of lattice artifacts. This gives a precise lattice definition of 
current and charge density independent of lattice size, and independent of the 
continuum limit. 

In the Abelian projection SU(2) link variables are parametrized by U(l) links 
and charged coset fields. The latter are normally discarded in Abelian pro- 
jection, as are the ghost fields arising from the gauge fixing procedure. Since 
the remainder of the SU(2) infrared physics must arise from these, an un- 
derstanding of their role is central to completing the picture of full SU(2) 
confinement. 

In the maximal Abelian gauge the supposedly unit charged Abelian Wilson 
loop has an upward renormalization of charge due to this current. A localised 
cloud of like polarity charge is induced in the vacuum in the vicinity of a 
source, producing an effect reminiscent of the antiscreening of charge in QCD. 
In other gauges studied, the analogous current is weaker, and acts to screen 
the source. 

We show that this current can be quantitatively written as a sum of terms 
from the coset and ghost fields. The contribution of the ghost fields in the 
maximally Abelian gauge in this context is found to be small. The effect of 
the the Gribov ambiguity on these currents is argued to be slight. 
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5.1 U(l) 



We first introduce and review the method due to Zach et al (|86|), in pure U(l) 
theories 

Consider a shift in the U(l) link angles in the partition function containing a 
Wilson loop source term 



Since the Haar measure is invariant under this shift, Zy/ is constant in these 
variables. Absorbing the shift into the integration variable and the taking the 
derivative 



Z w — 0, 



we get the relation 

= J [d0] {cos6 w - sm6 w /3A M sin 0^) e PE c ™ 9 ^. 
This can be cast into the form 



_ (sin^ sine,*,, 
\^^lw- {cQs6w) -eO x , Xw -J u . 



This is a well known technique to generate exact relations between Green's 
functions that is useful in generating Ward identities, or Schwinger-Dyson 
equations, or in this case we denote as Maxwell-Ehrenfest relations. We use 
the term Ehrenfest because it is the expectation value of what is normally a 
classical extremum of the path integral - an Euler Lagrange equation. 



5.2 SU (2) no gauge fixing 



Before addressing the full problem we first generalize from U(l) to SU(2) 
without the complication of gauge fixing. 



Z w ({U s }) = J[d(UU s )] W 3 (U) e 



PS{U) 
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The size of the source is irrelevant so we choose it to be the simplest case, i.e. 
a plaquette: 



W 3 ee -Tr{UWUUia 3 ). 
We choose the shift to be in the 3 direction 



Ufao) = (l - ^e M (x )<x 3 
The invariance 

d 



Z w = 



gives the Ehrenfest relation 



For 13 = 2.5, f3 (W 3 S^) = 0.0815(2), and (W) = 0.0818(1), and the difference 
= 0.0003(3). 

The notation denotes an e derivative(f45|). The denominator is just the 
ordinary plaquette. 



W ee ]-Tr{U ] U ] UU) 



To cast this into the form of Maxwell's equation we decompose the link into 
diagonal and off-diagonal parts 



U tx {x)=D lx {x)+O lx {x) 



Further we simplify notation 

(w 3 -..) 



)w — 



(W) 



We then group terms involving the diagonal part on the left and take all terms 
having at least one factor of the off-diagonal link to the right. 
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P ( s i*)w\ U=D - - ( A » F ^)w ■ 
Finally, note 

1 

e 

giving the final form of the Ehrenfest relation 



r 2_ jstatic 

°x,x w g ' 



(\F, u ) w = (j^) w + J^. 



This then tells us how to choose a lattice definition of field strength that 
satisfies an Ehrenfest relation: 



5.3 Gauge fixed SU (2), U(l) preserved 
The effect of gauge fixing gives 



Zw{{U s }) = J [d{UU s )\ W 3 (U) A FP 5[F] e^ u \ 
where we have introduced 

i = a fp [ nwn^ 1 ^)], 

j,y i,x 

and integrated out the g variables in the standard way. 

In this case Zy? is not invariant. The shift is inconsistent with the gauge con- 
dition. However, it is invariant under an infinitessimal shift together with an 
infinitessimal 'corrective' gauge transformation that restores the gauge fixing 



G(x) = (i - l -n(x) ■ <r) . 



Use of the invariance of the measure under combination of a shift and a 'cor- 
rective' gauge transformation we obtain 
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d d7] k (z) d 

de^zo) de t*( z o) drjk(z) 



'W 



In shorthand notation^, the Ehrenfest relation reads 



3)^ 



+ (W 3 ) M 



+ W 3 



(A 



FP)ii 



A 



FP 



(A 



FP)ti 



A 



FP 



+ /3S, 



0. 



(10) 



Gauge fixing has introduced three new terms: 



comes from the corrective gauge transformation acting on the source 



which is U(l) invariant but not SU(2) invariant. 



(Afp)p 



is due to the shift of the Faddeev-Popov determinant. 

is due to the corrective gauge transformation of the Faddeev-Popov 
determinant. 

The latter two derivatives are subtle. The key is to first consider the constraint 
up to first order in the shift and the corrective gauge transformations 



„. - dFAx) , , , ^ dFAx) , . , 

and then define the Faddeev-Popov matrix as a derivative of the corrected 
constraint. 



d 

drjAy) 



ijkhi i de ( z o 
Finally we evaluate the derivative using 



dFAx) , , . T-^dFAx) , . , 



£ = TriM-'M^ 



A check of the Ehrenfest theorem is given in Table 1. Some of the individual 
terms on the right hand side require a 2N x 2N matrix inversion, wherer N is 
the lattice volume. Hence to test the result numerically, we chose as small a 



4 SeeRef.f 
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lattice as possible. The result does not involve the size of the lattice which is 
4 4 in table 1. The last column employees a different source. The links making 
up the plaquette are replaced by the diagonal parts only. 

By separating the links into diagonal and off-diagonal parts we get the final 
form of the Ehrenfest-Maxwell relation. 



+ J, 



static 



+ (J, 



T FP 



+ u 



FP 



The first term in the current comes from the excitation of the charged coset 
fields, the static term has an extra non-local contribution coming from the 
corrective gauge transformation, and the last two contributions are from the 
ghost fields. 

These terms give a non vanishing charge density cloud around 
a static source. The lefthand side can be used as a lattice 
operator to measure the total charge density. 



5.4 Abelian point charge has cloud of like charge 



As a simple application we use this definition of flux to calculate div E on 
a source and the total flux away from the source. In Table 2, we see that 
the integrated flux on a plane between the charges plus the integrated flux 
on a back plane of the torus is larger than the div E on the source. The 
interpretation is that the bare charge is dressed with same polarity charge by 
the interactions and the neighborhood has a cloud of like charge. Hence there 
is ant i- screening. This charge density has contributions from all terms in the 
Ehrenfest relation. 



5. 5 Summary 



We have exploited symmetries of the lattice partition function to derive a set of 
exact, non-Abelian identities which define the Abelian field strength operator 
and a conserved electric current arising from the coset fields traditionally 
discarded in Abelian projection. The current has contributions from the action, 
the gauge fixing condition and the Faddeev-Popov operator. Numerical studies 
on small lattices verified the identity to within errors of a few per cent. We 
have found the Faddeev-Popov current in particular to be unusually sensitive 
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to systematic effects such as low numerical precision and poor random number 
generators, but the origin of any remaining, subtle biases, if they exist, is not 
clear; we have already considered all terms in the partition function. 

In a pure U(l) theory the static quark potential may be measured using Wilson 
loops that correspond to unit charges moving in closed loops, as demonstrated 
by \(A~F vl j)\ = 5w- In Abelian projected SU(2) the same measurements in 
the maximally Abelian gauge yield an asymptotic area law decay and a string 
tension that is only slightly less than the full non-Abelian value. In other 
gauges it is not clear that an area law exists — certainly it is more troublesome 
to identify. 

We have seen that in the context of the full theory the Abelian Wilson loop 
must be reinterpreted. The coset fields renormalize the charge of the loop as 
measured by \{A~F ut± )\ and charge is also induced in the surrounding vacuum. 
Full SU(2) has antiscreening/asymptotic freedom of colour charge, and in the 
maximally Abelian gauge alone have we seen analogous behavior, in that the 
source charge is increased and induces charge of like polarity in the neigh- 
bouring vacuum. Whether this renormalisation of charge can account for the 
reduction of the string tension upon Abelian projection in this gauge is not 
clear. In other gauges, where Abelian dominance of the string tension is not 
seen, the coset fields appear to have a qualitatively different behavior, acting 
to suppress and screen the source charge. 

In conclusion, the improved field strength expression defined by the Ehrenfest 
identity does not coincide with the lattice version of (|I7|) of 't Hooft's proposed 
field strength operator flUl) . The Abelian and monopole dominance of the 
string tension invites a dual superconductor hypothesis for confinement. If 
this is to be demonstrated quantitatively such as by verification of a (dual) 
London equation then a a careful understanding of the field strength operator 



is required. The Ehrenfest identities may provide this (|87j) 
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7 Figure captions 



Fig. 1. (figsu2.ps) (From ref. (|20|)) p vs. (3 for SU(2) gauge theory. The peak sig- 
nals deconfming phase transition, here monopoles are defined by the Abelian 
projection on Polyakov line. 

Fig. 2. (fig7.tex needs fig7 l.eps, fig7 2.eps and epsfig.sty) (From ref. (0)) 
Transverse profile of the electric field and curl of the monopole current in the 
mid plane between a static qq pair in the maximal Abelian gauge at finite 
temperature for a confining (left) and unconfining (right) phase. 

Fig 3. (curlE.ps) (From (]44])) Check of the dual Ampere law in the dual vortex 
profile. E is the electric field, k is the monopole current. 



Fig 4. (fit3.ps) (From (44])) Fit of the electric field vortex profile E and the 
tangential component of the magnetic monopole current kg to the Ginzburg 
Landau theory. 
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Source: 

Ehrenfest term 



W 3 



W 3 (U D) 



(W 3 )» 



0.65468(10) 0.63069(20) 



0.06095(7) 0.04463(4) 



0.00127(21) 0.00132(50) 

0.00529(3) 0.00564(3) 
-0.72246(68) -0.68275(50) 



Zero -0.00026(77) -0.00045(64) 

Table 1 

Terms in the Ehrenfest relation, Eqn. (terms). The column labeled 7 W 3 corresponds 
to the source described in the text. In the second column the links are replaced by 
the diagonal arts of the links in order to test a second example. The theorem gives 
zero for the sum. (3 = 2.5, 4 4 lattice 
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divE(cl.pt. charge) divE (on source) 



total flux 



10.0 

(almost 
classical) 



0.1 



0.1042(1) 



0.0910(8) (mid) 
0.0148(8) (back) 
0.1092(8) (total) 



2.4 0.4166 0.5385(19) 0.7455(70) (mid) 

0.0359(72) (back) 
0.7815(95) (total) 

Table 2 

divE normalized to ^ for a classical point charge. Source is a 3 Wilson loop. divE 
measured on a source. Electric Flux measured on the midplane centered on the 
Wilson loop. Also included is the flux through a plane on the far side of the torus, 
and the sum being the total flux. 8 4 lattice, 3x3 Wilson loop. 



28 



0.008 
0.007 
0.006 
0.005 
0.004 
0.003 
0.002 
0.001 




0.0 



■100.0 



-200.0 



-300.0 



-400.0 



1.0 



2.0 



3.0 
P 



4.0 



5.0 



0.03 
0.025 

0.02 
0.015 

0.01 
0.005 




0.40 



0.30 



0.20 



0.10 



0.00 



-0.10 







d=3a 




P=2.25 


o Electric field 


- ® 




\ -Curl J 




© 






A © 


® a 

* ® ffl ; 


- A 4 4 * 



0.0 1.0 



2.0 

r/a 



3.0 4.0 



0.10 



0.05 



0.00 



-0.05 



d=3a 



(3=2.40 



Electric field 
A -Curl J 



00f 



0.0 1.0 



2.0 

r/a 



3.0 4.0 
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